Theory of Magnetic Effects on the Noise in 
a Germanium Filament 

By HARRY SUHL 

(Manuscript received October 10, 1952) 

A magnetic field will influence the current noise in a germanium fila- 
ment. This fact bears out the hypothesis that at least part of the noise arises 
from minority carriers emitted in random bursts and recombining at the 
surfaces. A quantitative theory of this effect is given. 

INTRODUCTION 

Iii a series of fundamental experiments, H. C. Montgomery 1 has es- 
tablished that minority carriers play an important part in the current- 
noise associated with semiconductors. He found that on the one hand, 
the noise voltage is usually proportional to the biasing current, suggest- 
ing fluctuations in the conductivity, and hence the carrier concentration. 
On the other hand the spectrum of the noise suggested a rather coarse- 
grained time variation, not likely to be caused by fluctuations in the 
normal carrier density. One might conclude, therefore, that the noise 
is caused by a distribution of sources emitting or absorbing minority 
carriers in random bursts. Such carriers would be subject to the same 
laws of motion and of recombination as intentionally injected carriers. 
Montgomery was, hi fact, able to verify that the noise along a filament 
showed marked correlation over a distance roughly equal to that through 
which minority carriers could drift in the biasing field before recom- 
bination. 

W. Shockley has pointed out another corollary of this theory: A mag- 
netic field transverse to the filament should have a pronounced effect 
on the noise. This conclusion, too, Montgomery was able to verify 
experimentally. 1 His results are in good qualitative agreement with 
theory. Complete quantitative agreement was perhaps not to be ex- 
pected, since technical difficulties prevented attainment of the idealized 
conditions assumed by the theory. This paper gives an account of that 
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theory. On it are based the computed curves in Montgomery's paper 
showing the change in noise power with magnetic field. 

To see how such a change comes about, we imagine the magnetic 
field applied normal to one pair of the long faces of a rectangular fila- 
ment. This field, and the longitudinal drift current used to measure the 
noise, yield a sidewise thrust on the carriers, directed at right angles 
to the other pair of long sides. As a result the density distribution over 
the cross section is distorted, the minority carriers tend to accumulate 
near one of those sides, while the neighborhood of the opposite side is 
depleted. But for the usual conditions the recombination of carriers 
occurs mainly near the surface, and is proportional to their density 
there. Hence the magnetic field will change their lifetime. 2, 3 Clearly the 
amount of noise is dependent on the length of time carriers are able to 
contribute to the change in conductivity, that is, dependent on their 
lifetime. Therefore, the magnetic field should change the noise power. 
In simple extreme cases one can even make a semiquantitative argument 
for the maximum variation to be expected on the basis of such 
considerations . 

FORMULATION OF THE PROBLEM 

In order to make an exact calculation, we require a few preliminaries: 
The conductivity g is supposed to undergo a small time-dependent 
fluctuation Ag(t) about its mean value. 

The fluctuation arises from certain sources each of which, for macro- 
scopic purposes, may be considered to emit a noise-current J(t) of minor- 
ity carriers. Thus in a small time-interval dt' near t' the excess charge 
injected is J{?) dt'. This charge decays by recombination. Let r(t - t') de- 
note the fraction of carriers left over at time t(>t'). Then at time t 
there remains a charge r(t - t') J(t') dt' of the original injection. Now 
provided the excess density is small compared with the mean density, 
Ag(t) is proportional to the excess charge at time t, due to all the previous 
emissions added together. Therefore 

Ag{t) «: [ r(t- t')J{t')dt'. (1) 

J— 00 

In practice we do not literally plot Ag(t) as a function of t, but rather 
its frequency component Ag(J) in a narrow range df of frequencies near/. 
In other words, we single out for observation the contribution to Ag 
from that part J (J) of the injected current Jijt') which varies as e~ 2r,ft 
Suppose now that 1// is large compared with the time over which r(t) 



MAGNETIC EFFECTS ON NOISE IN A GERMANIUM FILAMENT 649 

is appreciably different from zero (that is, let 1// be much greater than 
the lifetime). Then, in the integral (i), r(t — t') will have gone from 
unity to zero long before J(f)e~ 2rift has changed appreciably from its 
value at t' = t. Therefore, for purposes of observation at frequencies 
much smaller than the reciprocal lifetime, we can rewrite (1) as 

Ag(t) cc J(t) J r(t - t') dt' 

(2)* 

= J(t) [ r(t) dt. 
Jo 

The integral in (2) can be interpreted as the average lifetime of car- 
riers. For, by definition, the rate of recombination at time t is — dr(t)/dt, 
so that — (dr/dt)dt is the number of carriers recombining between time 
t, t + dt. Hence the average lifetime is 

r = ~C l ^r dt = - WW" + C r{t) dt 

= J r(t) dt 
Jo 

since tr(t) —* as t — » oo . 

If 1// is not large compared with r one cannot simplify the integral 
(1) in this way. One then has to consider separately each frequency 
component Aglf)e~ 2 ' i/l due to the current J(f)e~ 2ri/t . Then 

*g(f)e- 2Ti " = J(f) f lit- t')e- 2Tif " dt' 

= J(f)e- 2rift C r{t!)e 2rift ' dt' 
Jo 

or 

Atftf) - J(f)r(f) 
where 

t(/) = I" r(t')e 2 ' if " dt'. 

Jo 

The calculation of t(/) is more complicated than that of t = r(0), and 

* From here on the equality sign will replace the proportionality sign. The 
resulting change of units is of no consequence in the final results which are only 
concerned with ratios of conductivity modulations. 
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at the present time the experimental situation does not call for refine- 
ments of this kind. Therefore we shall restrict ourselves to the calcula- 
tion of T. 

To evaluate t it is not necessary to consider a time-dependent case 
at all. In our experiment, it is the mean square conductivity fluctuation 
which is actually observed. Hence from (2) 

<Ag 2 > = <J\t)> r 2 . 

If the emission processes are stationary in time, <J"(t)> is time in- 
dependent: 

<Ag 2 > = <Jl> t. 

Now t can be written as 

[ r(t - dt', 

J— 00 

which is simply the total concentration at the present time t due to a 
constant injection from time — <*> to the present. 

Therefore the problem is reduced to finding the total carrier concen- 
tration in the filament due to a distribution of sources of constant 
strength y/ <Jl> . 

Let w(x, y, z; x x , iji , Zi) denote the carrier concentration at x, y, z 
due to a steady unit source at x\ , yi , %\. Then the total carrier con- 
centration is 

t(x-i , ?/i , *) = ] wfo y, z; Xi , yi , zO dx dy dz. 

The reason for the dependence on x\ , </i , Zi , is that the recombination 
process takes place largely on the surface. Therefore a source near the 
surface will yield a smaller concentration than one well inside the fila- 
ment. (Volume recombination will be neglected throughout this paper.) 
The mean square conductivity modulation due to many statistically 
independent sources at x r , y, , z T (r = 1, 2 • • • ) is then 

<(A0) 2 > = 2r'(a? r , Vr , *r) <J\Xr , Vr , Z r )> . 

The behavior of w is governed by the diffusion equation, subject to 
the boundary conditions expressing the recombination process, and sub- 
ject to a suitable singularity at x r , y r , z r , expressing the injection of a 
unit current. But in two and three dimensions the solution is not avail- 
able in closed form, or at any rate not in terms of the elementary trans- 
cendental functions. The infinite series for the solution is not easy to 
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handle computationally. It is therefore desirable to simplify the experi- 
mental conditions to a point where the problem becomes almost one- 
dimensional. A solution for w can then be found in closed form. 

Consider a very long uniform rectangular filament with one pair of 
sides very much wider than the other pair. Suppose that the y and z 
directions are respectively parallel to the wide sides and to the length 
of the filament (Fig. 1). 

Consider sources located anywhere on a plane x = £, which is parallel 
to the wide sides of the filament. If the recombination properties of the 
filament are uniform in the y — z directions, the lifetime due to a unit 
source anywhere in that plane is independent of the location of the 
source on that plane, and depends only on £. Hence the conductivity 
modulation due to sources of strength 



J(£, Vt , Zr) 



(r= 1,2-..) 



in that plane is simply 



T(S)T,J&yr,Zr) 



. 
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Fig. 1 — Geometry of the filament, and disposition of the fields. 
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and is the same as that due to an infinite source of strength 

22 J(t, Vr , Z r ) 

T 

uniformly distributed over the plane x-£. 

But the density w due to such a source will be a function of x and£ 
only, and will be the solution of the one-dimensional diffusion equation. 
Hence for a geometry approaching that of figure 1 sufficiently closely, 
the problem is one-dimensional. 

The Evaluation of r 

We now have to write down the one-dimensional diffusion equation 
in the presence of a magnetic field along Oj, , which combines with the 
drift velocity of the carriers so as to force them towards one of the 
surfaces x = ± a. If F x is the effective field arising in this manner, and 
D is the diffusion coefficient, the equation is 

which expresses the fact that the diffusion current — Dq dw/dx plus the 
drift current qnF x w must be constant since the carrier density cannot 
build up indefinitely, n is the mobility of the minority carriers: /m for 
electrons, n P for holes. As is shown elsewhere 2 the effective field F x is 
given by 

F x = (fi n + e p ) E z = 0E, , 

where E z is the biasing field causing the drift current, 0„ , 6 P are the Hall 
angles for electrons and holes, respectively. If y. n , n P are the electronic 
and hole-mobilities, and if the magnetic field is not too large, 

e = e n + e p = lcrVn + v v )H, 

where H is in oersteds, and the mobilities are in cm 2 /volt-second, and 6 
is in radians. (Strictly speaking, the diffusion current is not in the direc- 
tion of the density gradient when a magnetic field is present. As the 
result mixed derivatives 

a w 



dx dy 



occur in the diffusion equation. But in the reduction to one dimension 
these terms integrate out. All that remains is a small correction to D, 



MAGNETIC EFFECTS ON NOISE IN A GERMANIUM FILAMENT 653 

negligible for ordinary values of H.) It is convenient to specify a dimen- 
sionless parameter in the same notation as H. C. Montgomery. 

2a l iF x = 2aBE z 
D ~ D/n 

By the Einstein "Relation D/y. = kT/q this may be written 

_ 2aBE z 

" kT/q 
where q is the absolute value of the electronic charge. $ is the ratio of 
the voltage corresponding to the transverse field to the thermal voltage 
kT/q. In terms of <1>, equation (3) can be rewritten 

d"w _ «I> dw _ _ , v 

rtV ~ 2a tfr ~" w 

The integral of this equation has the form 

w = Ae* (l/2a) + B (5) 

where ^4 and B are two constants. Because of the existence of a singu- 
larity at x = .r , say, the constants A, B take on different values for 
x <x and x>x ■ To see what these values are, we first write the solu- 
tion (5) in the form 

Wl m ^ ie *<*-*o>/2a + ^ x > x ^ 

W2 = A 2 e* lx ~ X0)l2a + B 2 x <x . 

At x = x the density w must be continuous. Hence 

Ai + Bi = Ai + B 2 (6) 

Further, the discontinuity at x must be such that the difference of the 
currents on the two sides is just unity, the strength of the injected 
current. Now the total current is 



/dw 4> \ 
— D [ -, — w ) 

\dx 2a ) 



(i.e., the diffusion current plus the drift current), and w is continuous 
at .T . Hence the difference between the current on the two sides of x Q is 



l--DLun(fe) -(^. 

A-o \\dx / Xo +h \ dx / Xo -h 

= I- D(A 2 - A x ). 
2a 



(7) 
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So far we have two relations (6) and (7) between the four constants 
Ai, A if Bi, Bi. To determine these constants, we need two more 
equations. These are supplied by the boundary conditions. We assume 
that the recombination rate is proportional to the density at the bound- 
aries (x = ±a). The recombination then has the formal appearance of 
a current through the surface, the factor of proportionality s playing 
the part of a "recombination velocity." 2, 3 That current through the 
surface must equal the current arriving at the surface from the interior 
of the filament, under steady-state conditions. The boundary condi- 
tions are thus: 

- D fe-l w 'L. =s,u,(+o) ' 

-D[-P -77-w*) = -SiW(-a), 
\ ax 2a /x—a 

where S\. , S2 are the surface recombination velocities of the faces x = 
+a, — a respectively . The minus sign on the right-hand side of the second 
of these equations is due to the fact that the recombination current at 
x = — a is along the — x direction. Defining the recombination param- 
eters 

*. 2 - -g- 
and substituting the solutions W\ , w% in the last two equations we obtain 

B * = fcUrf*' 1 -"""" + B.) 

2 (8) 

Equations (6), (7) and (8) suffice to determine A\ , At , B% , B% . We 
easily find that 

, -*(W)/i 2a ai (l - a,e* ll+XoM ' 2 ) 

Al " ~ e m A 



o„ ft I — *(l-x /a)/2\ 

, , +$(l+z /a) /2 *a «2U + aie ) 

Ai = + e m a 

2a (1 - a2 6* (1+ * 0/ ' )/2 ) 

Bi ~ ~£5 a 

2a (1 + a ie -* (, - Io/o>/2 ) 
B * = ~*D A 



(9) 
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MAGNETIC FIELD PARAMETER, $ 

Fig. 2 — Noise change for a filament with equal surfaces as a function of * for 
various ^. Surface generation. 



where 



«i - ( 2 " ft W> I a * = \2 +^ ! W 2 ; 



_ -, -*(l-xo/<»)/2 



A = aie 



+ «2e 



*(l+«o/o)/2 



Now we are in a position to calculate t b (x ), the life-time in a magnetic 
field # (contained in $>), due to an emission at x Q : 

/x /.+o 

w%{x, x ) dx + / u-ifoxo) da: 
a •'in 



= (£i + 5 2 )a + (£, - BJxn + ^ (A, -Ai) 

2a r^ g *(l-x /a)/2 _ ^ 2g -*(l+io/<»>/2j 

Special Cases 

With the help of (9) and (10) we can now examine special cases. 

1. Surface Emission 



10 
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When only the surface x = +a is emitting, we set x = -f-a in (10) 
and get 



r H {+a) = 2B 2 a + ^ A 2 (l - e~*) 



4£ 



(1 + «i) 



l — « 2 — (1 — C ) 



2a 



a« 






$ 



ai + a 2 e 



Similarly, when only —a is emitting, we get 



2 «1 



e — 1 



T/, 



("«) = " ^7" 



<I> 



- 1 



D\p> ot\e * + (X2 



But in an actual experiment, both faces will be emitting, with mean 
square strengths (J\), {A) say. The quantity that is then measured 




-20 -16 



3-404 8 

MAGNETIC FIELD PARAMETER, $ 

Fig. 3 — Contribution to the noise change from a unit source at the plane 

x = +a. 
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is the ratio 



N H = 



Aflri 



(Jl)rl(+a) + (Jl)rl(-a) 



AgU (Jl)rU (+«) + <^>ri= (-a) ' 



:n) 



Therefore we also need the lifetimes at zero H (that is, zero $)• But at 
$ = 0, the r's are indeterminate, and we therefore have to take limits. 
Expansion in powers of $ shows that 



It r B (+a) = £- 

s— o JJyi 



1 + ^ 

V2 



1 + SV^! + ^ 2 



It T H (-a) = 
*-«o 



Thus we finally get 






-f 1 

2\*i 



a 



1 + 



1 



1+14- + 4 



*1 2 




2 


aie * -\-a 2 


2 



(12) 



<£[> 
« 



'+£ 



2 , <^> 



1 + i- 



There remains one small difficulty in the way of comparing experiment 
with theory: We do not know (Jt), (Jt)- As suggested by H. C. Mont- 
gomery, we are able to overcome this difficulty as follows : We first draw 
a number of curves of 



rl(+a) 



ri(-a) 



rl(+a) ' rl(-a) 



versus $, 



for various sets of parameters (^1 , ^2.) (See figures 3, 4). Then we 
contrive to match a linear superposition 



C\ 



Tfl(+a) , . Ta(-a) 



To 



(+a) 



+ C2 



To 



(-a) 



where Ci + d = 1 to the experimental curve. This will be possible only 
for one particular set of values (^1 , \f/ 2 ). From Ci and Ci = 1 — Ci and 



658 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1953 




-8-4 4 8 

MAGNETIC FIELD PARAMETER, $ 



-20 -16 -\i 

Fig. 4 — Contribution to the noise change from a unit source at x = —a. 



from (^i , fa) we can determine the ratio of (Jl) (Jl), for, having matched 
the experimental curve, we know that 

rg( + g) 4. /1 — ) T *(~" g ) - (^i) T g(+ q ) + (Jl)r%( — a) 
Cl rl(+a) + Cl) rj(-a) (J?>r5(+a) + <^>r!(-a) ' 

which is satisfied for all $ if 

ft ^ <jj> 

r^(+a) <JOTo(+a) + (Jl>'o(-a)' 

(1 ~ ci) _ (jl) 

rj(-a) < J?)rJ(+a) + < JiJfR-a) " 

These equations are self-consistent and give 
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ft _ Tl(+a)(Jl) __ (Jl) \ l + TJ $ 



1 " ft rl(-a)(Jl) (Jl) / ~ JV tf ' 

Thus the match of theory to experiment actually provides information 
about the relative emissivities of the surfaces. 

When the surfaces are equal, yf/ x = \f/ 2 = $; (Jl) = (Jl) and the 
result 12 simplifies to 

.2 



_('+{)'+ 8 -■*# 



N "- r , # x. ' w 

1 + 2? COth 2, 

10 log A 7 * for this case is plotted in figure 2. 

It is interesting to see how N B for equal surfaces varies for small 
values of $. 2/$ — coth <£>/2 is regular at $ = and varies as — <I>/6 
there. Hence the initial variation is as 



AT 



Hence the curve rises or falls initially according as 

The noise therefore increases initially if 

5" 1 1 

yx — = > j or ^ > 6.9 approximately 

and falls initially if 

i> < 6.9. 

2. Volume Generation. 

At first sight it may seem that if volume generation is considered, so 
should volume recombination (detailed balancing). This is not neces- 
sarily so, since we are not dealing with an equilibrium situation here. 






060 
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Therefore the possibility of volume generation and surface recombina- 
tion cannot be discarded. 

This case is somewhat more difficult. Assuming all sources to be un- 
correlated and uniformly distributed throughout the interval (-a, +a), 
we have to square expression (10) and integrate it from x = —a to x = 
+a in order to find (Agl). (We suppose that all the sources are of equal 
strength (J 2 )). Substituting the values of the A and B from 9 and 10, 
we get, after some obvious cancellations 






where 



S = 2a 



(i + »(&!#) 



aje 



,-*/2 



+ «2e 



•/a 






, -*/2 +*/2 

2 a ie -* /2 + a 2 e+* /2 



aiota sinh 



$ 



- 2 



Hence 



4a 
4> 2 Z> 2 



me'*' 2 -f- a 2 e* /2 
t 2 ,(x ) dx 

a 



- 1 



? a 3 + 23 .V + ^sinh| 
3 <$> 2 



(14) 






Before proceeding with this general case, we first consider the limiting 
case i^i = fa = «>, when ai = — 1, a 2 = +1. Then 



T = -a coth 



* 



5 = 



sinh 



$ 



and 



f +a a 8a 

L Th{Xo} ~ m 



&D 2 



1 *•• r C ° th l,8 

3 + coth 2 ~ G -J - + * 2 
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To find N H = (Ag 2 n )/(Ag 2 H=0 ) we need the limit of ft a a T 2 H as $ -> 0. 
After some tedious algebra, we find this limit to be 



Tff=0 — 



4a 
15Z) 2 



so that 



\Nn)t i.*2=°° — 



30 

$2 



1+ nth'* GCOth 2^8- 
3 + C ° th 2--^-+^ 



(15) 



In the general case we can again take the limit of (14) as $ — ► in 
order to determine N H , but this would be too tedious. Instead, we 
solve the diffusion equation directly when $ = 0. The equation is then 
simply 

s- 

and the solution subject to the correct boundary conditions and allowing 
for a steady unit injection at .r is 



u>i = Ai(x - .T ) + B 
Wi = A 2 (x — x ) + B 



x > x 

X < .T . 



where 



^i - -is* 



D^ 



1 I (1+j 



/ OAl + ^2 + 2M$> 

(h + fo + 2^ 2 ), 
B = I [l + * (l - 5)] [i + ^ 2 (i + ]*)]/(* + ^ + a**). 



A.-^ri+^ri-j 



We now have 



w dx = / w 2 dx + / i^ dx 

a J- a Jx 

* + *] + 2aB. 



= -- | 2ax (A l + A 2 ) 4 



From this we can compute 



/• +a 2 

/ ri_o(a: ) ^/x , 

J— a 
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16 



Fig. 5 — Noise change in a filament with uniform volume generation and equal 
surface recombination rates. 

but the result is still rather complicated unless ^i = ^2 = $• If we there- 
fore restrict ourselves to that case we find 



Ax + A 2 = - 



x 



and 



and 



D(l + $) a' 

2Z) L * 1 + * ^ J 



Dt b -o(x ) = a ( - + 7 ) - 77 



r*-w— s<i+6+j)> 
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(Note that as \f/ — * «> , this quantity tends to 4a 5 /15Z) 2 , as before in our 
limiting process.) 

When \j/i = ^ 2 = ^, S and T also simplify somewhat, and the result is 






= 2 



where now 



3 + 2T ~ + ¥ U ~ T J smh 2 ~ ¥ cosh 2 + T smh * 



**kh)+2 



1+ ] 

5- * 



cosh = + sinh 5 J 



$ -I' 

— + coth - 

r . i , 1 ■ i V 2 



-(•♦& 



'l 1 <1> 

+ 2* COth 2 



An alternative form is 

4 



Ne = 



45 



i+G+,7±+A0+})6-*i) 

-, (16) 




v, 2 coth I 

P 4- coth?) 4-- - 



where 



2<£ V V $ 



«¥> <J> 



a result which correctly tends to (15) as $ — » 00 . 10 log N B for various 
\f/ is shown in Fig. 5. 
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